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Abstract - In the present study, we investigate patterns in the postural sway that character- 
ize the static balance in human beings. To measure the postural sway, sixteen healthy young 
subjects performed quiet stance tasks providing the center-of-pressure (COP) trajectories. From 
these trajectories, we obtained the COP velocities. We verified that the velocity distributions ex- 
hibit non-normal behavior and can be approximated by generalized Gaussians with fat tails. We 
also discuss possible implications of modeling COP velocity by using generalized Fokker-Planck 
equations related to Tsallis statistics and Richardson anomalous diffusion. 
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In daily life, we perform several tasks that demand 
a constant maintenance and control of postural balance 
(e.g., walking, dancing, carrying shopping bags, and so 
on). From a biological point of view, the quiet up- 
right stance requires of the body the control of the mus- 
culoskeletal system. Such control is possible when the 
central nervous system processes visual, vestibular and so- 
matosensory information, and uses it to select and control 
the appropriate postural response [1-5]. 

The study of postural sway by using center-of-pressure 
(COP)-based measures has been useful to characterize 
the dynamics of the postural control system associated 
with maintaining balance during quiet standing. The dy- 
namics of COP position, which is frequently related to 
a random walk, has been investigated using a variety of 
techniques [6-16]. In particular, some studies have focused 
COP velocities instead of trajectories. It has been argued 
that COP velocity is most accurate than position or ac- 
celeration to modulate human quiet stance [17-19]. Sev- 
eral characteristics of COP velocity, including temporal 
correlations, have been explored [19-25]. However, inves- 
tigations related to velocity distributions are scarce. For 
instance, there are statistical significance tests indicating 
that COP velocity distributions are non-Gaussian [26,27], 
but more details about the shape of these distributions 
have not been studied. The present study is mainly di- 
rected to investigate COP velocity distributions. 



In this work, sixteen healthy subjects participated in the 
experiment (eight females and eight males between twenty 
and twenty-eight years old) . After the ethics committee of 
the Universidade Estadual de Maringa approved the study, 
all participants gave their informed consent prior to their 
participation. All participants received oral and written 
instructions about the experimental protocol. They signed 
a commitment consenting before performing the tasks of 
the experiment. 

The volunteers stood barefoot and were asked to 
maintain quiet upright posture on the force platform at 
a distance of about 1.5 meters from a fixed point in the 
horizontal plane with their legs oriented in parallel at hip 
width apart. Each subject was submitted to this task 
for 60 s, and the same was repeated 10 times in an alter- 
nating fashion, i.e., participants took turns among them- 
selves in order to alleviate signs of fatigue. Thus, a total 
of N = 160 time series for this experiment was obtained. 
Data were collected using a force platform (EMG SYS- 
TEM do Brasil) to record the ground reaction forces with 
an acquisition rate of 100 Hz. 

The collected COP trajectory time series were filtered 
by a low-pass filter with a cut-off frequency of 20 Hz. This 
procedure eliminates high-frequency fluctuations. We 
then considered a sampling frequency of 50 Hz. The ve- 
locity series were obtained by differentiating the position 
series. The velocity time series were visually examined and 
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Fig. 1: (A) A typical 60s center-of-pressure (COP) trajectory, where X t and Y t correspond to the mediolateral and anteroposte- 
rior directions, respectively. (B) The time series of the COP position filtered by a low-pass filter (top), Yt, and its corresponding 
velocity series (bottom), AY*. 
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Fig. 2: Statistical aspects of COP velocity series. Values of (A) standard deviation, (B) skewness and (C) kurtosis of the COP 
velocity for each subject (anteroposterior direction). The dotted lines represent the expected value for a normal distribution. 
(D) The moments of the standardized COP velocity distribution (black) in comparison with the Gaussian moments (gray). 



it was decided that the first 1.5 second of each series would 
be removed to avoid anomalies related to the destabiliza- 
tion of the subject at the beginning of the experiment. 

A representative COP trajectory is displayed in fig. 1(A) 
for a subject. Although COP presents a two-dimensional 
behavior, we chose to emphasize the investigation on 
its velocity in the anteroposterior direction (forward and 
backward body movement) as each individual has a higher 
propensity to imbalance. By using Y t and Y t +i to repre- 
sent two consecutive terms of the COP position series in 
the anteroposterior direction, one has AY t = Y t +i — Y t as 
an element of the series of successive increments. The in- 
crement AY t will be referred to as COP velocity (Alt 
is proportional to the velocity). Figure 1(B) shows a 
typical position series and its corresponding velocity se- 
ries for the anteroposterior direction. Moreover, stan- 
dardized increments (standardized velocity) are written 



as Ayt = (AYt — //)/cr, where \i is the mean value and a is 
the standard deviation of AY t . Points above six standard 
deviations were ignored in all series because it was under- 
stood that these data include involuntary stimulation of 
the subjects. 

The asymmetrical character and "peakedness" for each 
COP velocity series were studied. It was also adopted 
another way to quantify the probability distribution of 
our total data set which is given by the moments of the 
distribution. Furthermore, in the direction of suggesting 
candidates to model COP velocity distributions, two fam- 
ilies of non-Gaussians with fat tails related to generalized 
Fokker-Planck equations were considered. 

To obtain the following results, the procedure taken was 
to group the ten velocity series of each subject in a single 
data set. The standard deviation, a, of AY t for each sub- 
ject depends a lot of the subject characteristics (fig. 2(A)), 
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a result consistent with the literature [28] . The asymmetry 
measure of all COP velocity distributions does not present 
considerable deviations when compared to a symmetrical 
distribution as the normal one, z.e., each skewness is ap- 
proximately equal to 0 (fig. 2(B)). The skewness measure 
employed here is {(Ay t ) 3 ), where (• • • ) denotes the average 
over the standardized COP velocities of a subject. Nev- 
ertheless, all measures of the "peakedness" of the proba- 
bility distribution indicate a non- Gaussian behavior since 
the values of the kurtosis are significantly larger than 3 
(fig. 2(C)). Here, the kurtosis measure is given by {(Ay t ) 4 ). 

In order to overcome the variability of the standard de- 
viations, as in the cases of skewness and kurtosis analysis, 
and to try to identify more patterns related to COP veloc- 
ities, the following analysis is based on the standardized 
COP velocity Ay t . In this direction, consider the mo- 
ments jik of the distribution of standardized COP veloc- 
ity Ay tl with /ife = (\Ay t \ k ), where all the series were used 
as a single normalized data set. Figure 2(D) shows that 
the moments of the velocity distribution deviate from the 
Gaussian case as k increases. This result extends the kur- 
tosis analysis for the global data because is the kurtosis. 

In addition to the above analysis, tests of statistical sig- 
nificance were performed for data of each subject. The 
null hypothesis that the data is distributed according to 
the Gaussian distribution is rejected at the 5 percent level 
based on the Anderson-Darling and Cramer-von Mises 
tests. This analysis is consistent with previous studies 
where it was also indicated that COP velocity distribu- 
tions are significantly different from Gaussians [26,27]. 

Next, the COP velocity distribution is explicitly calcu- 
lated, i.e., P(Ay t ) is obtained by using all N = 160 series 
as a single normalized data set. From fig. 3, one veri- 
fies that P(Ay t ) remembers the normal distribution for 
\Ayt\ ^ 3. In contrast, for \Ay t \ > 3, a deviation from 
a Gaussian is identified. A direct practical implication of 
this finding is that large velocities occur with larger prob- 
ability than expected for a Gaussian process. 

Since COP movements can be viewed as random walks, 
candidates for COP velocity distributions could be re- 
lated to Fokker-Planck equations. In the usual case of 
one-dimensional free-particle motion, the Fokker-Planck 
equation for the velocity v is given by [29] 



dP d f^dP\ d , m 



(1) 



where the noise strength D is constant, the friction co- 
efficient 7 is also constant and P(v,t) is the probability 
velocity distribution at time t. In the stationary case 
(dP/dt = 0), with P(v) — > 0 for \v\ — > oo, this equation 
reduces to 

dP 

L>--+7^P = 0. (2) 
dv 

The solution of this last equation is the Maxwell velocity 
distribution, i.e., the Gaussian P(y) = N' exp(— jv 2 /2D), 
with N' being the normalization constant. As the 




Fig. 3: (Colour on-line) Standardized anteroposterior COP 
global velocity distribution, P(Ay t ). The dotted line rep- 
resents a Gaussian with zero mean and unit standard devi- 
ation. The continuous lines refer to a g-Gaussian (eq. (3)) 
with q = 1.244 (red) and a stretched Gaussian (eq. (4)) where 
s — 1.301 (blue), values obtained by using the maximum like- 
lihood estimation. 



COP velocity distribution exhibits non- Gaussian behav- 
ior, candidates for COP velocity distributions (generalized 
Maxwell velocity distributions) could be investigated from 
stationary solutions of generalized Fokker-Planck equa- 
tions. Two cases of generalized Fokker-Planck equations 
that lead to a significant improvement in the description 
of the COP velocity distribution, when compared with a 
Gaussian, are focused here. 

As a first generalization, consider D = D q P 1 ~ q with 
D q constant. This leads to a non-linear Fokker-Planck 
equation [30,31] which recovers the usual one when q = 1. 
The corresponding solution of eq. (2), when replacing v by 
Ayt, is 



P q {Ay t ) = N q [l-{l-q)bAy, 



(3) 



where N q is a normalization constant and b is a function 
of D q and 7. Note that the Gaussian distribution cor- 
responds to P q (Ay t ) in the limit q —> 1 and P q (Ay t ) has 
fat tails because P q (Ay t ) oc \Ay t \ 2 ^ l ~ q ^ for large \Ay t \ 
and q > 1. In addition, P q (Ay t ) w N q (l - b Ay 2 ) for 
Ay 2 <C l/b. Therefore, eq. (3) has two aspects of COP 
velocity: a central part resembling a Gaussian and long 
tails when q > 1. 

Equation (3) has been employed in many contexts, from 
physics and engineering to economics and biology [32,33]. 
For instance, it was applied with success to model the 
non- Gaussian velocity distribution of Hydra viridissima 
cells in two kinds of cellular aggregates [34]. Because 
the connection of eq. (3) with the Tsallis generalized sta- 
tistical mechanics, it is usually referred as g-Gaussian. 
This statistics has been a useful theoretical tool to dis- 
cuss anomalous behaviors, in particular, when related to 
non-Gaussians [32,33]. For instance, the usual Maxwell 
distribution is replaced by eq. (3). In this scenario, many 
aspects related to the usual (Boltzmann-Gibbs) statistical 
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mechanics have been generalized. An example is the use 
of the non-linear Fokker-Planck equations [30,31,35,36] re- 
placing the usual ones. 

As well as the non-linear Fokker-Planck equation was 
first introduced in connection with a non-linear anomalous 
diffusion equation (porous medium equation) [30] , another 
generalized Fokker-Planck equation can be related to other 
anomalous diffusion. In this direction, the same functional 
form of Richardson anomalous diffusion equation [37] is 
employed here but replacing the distance r by the veloc- 
ity v. This proposal corresponds to using D = D s \v\ s ~ 2 
in eq. (1), with D s and s being constants. In this case, 
the solution of eq. (2), a generalized Maxwell distribution, 
becomes 

P s (Ay t )=N s exp(-a\Ay t \ s ), (4) 

where v was replaced by Ay tl N s is a normalization con- 
stant and a is a function of s, D s and 7. Equation (4) 
reduces to a Gaussian for s = 2 and becomes stretched 
for s < 2. For this reason, eq. (4) is referred to here 
as stretched Gaussian. Analogously to the g-Gaussian, 
the stretched Gaussian is common in many scenarios, 
such as diffusion in fractals [38], two-dimensional turbu- 
lence [39], predator-prey encounters [40], atom deposition 
into a porous substrate [41], distribution of sound ampli- 
tude in music [42], and ergodicity breaking in heteroge- 
neous diffusion processes [43,44]. 

Since the standard deviation of Ay t is equal to one, 
standardized q- Gaussian and stretched Gaussian distri- 
butions must be employed. As consequence, one has 
b = 1/(5 — 3g), with q < 5/3, in eq. (3) as well as 
a = [sT(l + l/s)]- s / 2 r(3/s) s / 2 in eq. (4), with T(x) be- 
ing the Gamma function. By using the maximum likeli- 
hood method, the COP velocity distribution for our global 
data was approximated by a q- Gaussian and a stretched 
Gaussian. Figure 3 suggests that the use of g-Gaussians 
and stretched Gaussians improves substantially the ad- 
justment of the COP velocities in comparison with a 
Gaussian. 

The same procedure was repeated to analyze the ve- 
locity distribution for each subject, ie., P( % \Ay t ), with 
i = 1,2,..., 16. In this case, as previously pointed out, 
the ten series for each volunteer were grouped into a se- 
ries and each P^(Ay t ) was approximated by standard- 
ized g-Gaussian and stretched Gaussian distributions. An 
illustration of an individual COP velocity distribution for 
a representative subject is presented in fig. 4(A). The 
results of the individual analysis are shown in fig. 4(B) 
(1.165 < q < 1.309 and 1.179 < s < 1.488), indicating 
that the q parameter remains larger than 1 and the s pa- 
rameter is smaller than 2. Moreover, these parameters 
are close to their global values obtained in the previous 
analysis, q = 1.244 and s = 1.301. Tests of statistical sig- 
nificance (Anderson-Darling and Cramer-von Mises tests) 
indicated that more than 70% (80%) of 160 anteroposte- 
rior COP velocity series are not rejected at 1% level for 
g-Gaussians (stretched Gaussians). 
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Fig. 4: (Colour on-line) (A) Standardized anteroposterior COP 
velocity distribution, P(Ay t ), for a representative subject. The 
dotted line represents a Gaussian with zero mean and unit 
standard deviation. The continuous lines refer to a ^-Gaussian 
with q = 1.296 (red) and a stretched Gaussian where s = 1.196 
(blue). (B) Values of the q (top) and s (bottom) parameters 
obtained for each subject. The solid lines are distances between 
the experimental values and those obtained from a Gaussian. 

A similar analysis was performed to the mediolateral 
direction (left-right movement as X t in fig. 1(A)). One 
found q — 1.177 and s = 1.443 for the global velocity 
distribution. For the individual cases, 1.104 < q < 1.226 
and 1.296 < s < 1.619 were obtained. The same tests 
of statistical significance (Anderson-Darling and Cramer- 
von Mises) pointed out that more than 50% (60%) of 160 
velocity series are not rejected at 1% level for q- Gaussians 
(stretched Gaussians). 

In addition to the analysis of COP velocity distri- 
butions, we present below results concerning COP po- 
sition distributions. Following a similar procedure to 
the velocity series, the fitting parameters obtained for 
the global position distribution were q close to 1 when 
employing q- Gaussians and s close to 2 for stretched 
Gaussians for both directions, anteroposterior (AP) and 
mediolateral (ML) ones. Explicitly, these values were 
q = 1.01 (q = 1.059) and s = 2.023 (s = 1.928) for 
AP (ML). Skewness in both directions was near to zero, 
-0.02 (-0.07), and kurtosis 3.02 (3.22) for AP (ML) 
direction. For each individual, the values of g, 5, skewness 
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and kurtosis remained close to the global values. For 
the tests of statistical significance (Anderson-Darling and 
Cramer-von Mises) less than 18% (22%) of 160 posi- 
tion series were not rejected at 1% level for g- Gaussians 
(stretched Gaussians) in both directions. These results 
show that the non-Gaussian behavior for the position 
series is subtle or almost imperceptible. 

To conclude our discussion about COP velocity distribu- 
tion, the main achievements are summarized and possible 
implications of our findings are raised. Consistently with 
the literature [28], fig. 2(A) indicated that standard devi- 
ations exhibit large variability from a subject to another 
and therefore cannot be used as a typical characteristic 
of all subjects. To investigate more deeply aspects unre- 
lated to the standard deviation, a standardized COP ve- 
locity was employed here. In contrast to the investigation 
of standard deviations, the skewness analysis (fig. 2(B)) 
pointed out in the direction that COP velocity distri- 
butions are well approximated by symmetric ones. In 
addition to this common aspect for all subjects of our 
analysis, the kurtosis values indicated a robust fat- tail 
pattern for COP velocity distributions when compared 
with a Gaussian (fig. 2(C)) and reinforced by a moment 
study (fig. 2(D)). The non-Gaussian profile was also ver- 
ified by tests to decide if the velocity data comes from 
a Gaussian distribution, consistently with previous stud- 
ies [26,27]. In addition, for standardized velocities, very 
similar distributions were identified for all subjects. And, 
considering generalized Fokker-Planck equations for COP 
velocities, the non-Gaussians distributions with fat tails 
can be interpreted as generalized Maxwell velocity distri- 
butions. These distributions were well approximated by 
g-Gaussians and by stretched Gaussians (figs. 3 and 4). 

The empirical results obtained here for COP veloc- 
ity distributions, together with those about correla- 
tions [19-21,24], compose a useful guide to model COP 
dynamics. However, others studies are still required to 
obtain a broader scenario for COP velocity distributions. 
Future works could investigate the sensitivity of COP ve- 
locity distributions to changes in postural steadiness re- 
lated to age, neurological and orthopedic diseases or other 
pathologies affecting the postural control system. 

Concerning the models employed here, they represent a 
first view of the scenario of COP velocities. In the model 
related to g-Gaussians, the noise strength D is propor- 
tional to a negative power of the probability distribution, 
P 1 ~ q with q > 1, leading a non-linear stochastic process. 
For the case of stretched Gaussians, the noise strength 
is proportional to \v\ s ~ 2 with s < 2, that is, a heteroge- 
neous model. If one is interested in a model that accom- 
modates g-Gaussians and stretched Gaussians in a unified 
way [45], D oc P 1 ~ q \v\ s ~ 2 could be used. Another possibil- 
ity would be to consider other models that have stretched 
Gaussians as solutions. This is the case of the fractional 
Fokker-Planck equation [46,47], which is equivalent to con- 
tinuous time random walk models with scale-free waiting 
times. Such a model presents a degree of memory and 



could be useful in connection with some regulatory mech- 
anism in the maintenance of the upright standing. In turn, 
this Fokker-Planck equation and heterogeneous models ex- 
hibit weak non-ergodic behavior [43,44,48]. Motivated by 
a possible connection between these models and COP dy- 
namics, the time dependence of the mean square displace- 
ment, measured from single COP trajectories, could be 
investigated. This procedure would give direct experi- 
mental verification of an anomalous stochastic process and 
of a non-ergodic character. Furthermore, to employ this 
procedure is useful because probability distributions could 
not be the best way to identify anomalous stochastic pro- 
cesses. All these aspects illustrate possible perspectives 
for further analysis involving upright quiet stance data. 
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